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Abstract 

We present an analytical method for finding all the electronic eigenfunctions 
and eigenvalues of the aperiodic Thue-Morse lattice. We prove that this sys- 
tem supports only extended electronic states which is a very unusual behavior 
for this class of systems, and so far as we know, this is the only example of 
a quasiperiodic or aperiodic system in which critical or localized states are 
totally absent in the spectrum. Interestingly we observe that the symmetry of 
the lattice leads to the existence of degenerate eigenstates and all the eigen- 
values excepting the four global band edges are doubly degenerate. We show 
exactly that the Landauer resistivity is zero for all the degenerate eigenvalues 
and it scales as ~ (L = system size) at the global band edges. We also 

find that the localization length ^ is always greater than the system size. 
PACS numbers: 71.23.Ft,61.44.-n 
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In recent years, intense theoretical activities have been observed in the hterature towards 
the understanding of the electronic properties of quasiperiodic and aperiodic systems. It has 
been seen that apart from the usual generic features like the critical wave functions, singular 
continuous energy spectrum, scaling properties of the density of states, power-law behavior of 
electrical conductivity etc. these systems can also sustain many features characteristic 

of periodic as well as disordered systems. As for example, in a very recent work Macia and 
Dominguez-Adame [§] have reported a new class of energy eigenstates in the spectrum of ID 
quasiperiodic Fibonacci lattice having properties characteristic of extended electronic states. 
But Fibonacci lattice is the simplest and the most well-studied quasiperiodic model system 
since the discovery of icosahedral symmetry in Al-Mn alloy by Shechtman et al. |^ and the 
general consensus about its electronic structure is that all the eigenstates are critical with 
Cantor-set energy spectrum. This kind of diversity of the electronic and other properties 
are the general characteristics of almost every quasiperiodic or aperiodic system. A few such 
instances are the period-doubling 0,0], copper- mean |0, Rudin-Shapiro and Thue- Morse 
T^,[TT]] etc. lattices and in these systems critical as well as extended wave functions can 



coexist. Moreover, in some of these cases there is evidence of localized |T2| or chaotic wave 



functions ||T^. There are also notions like the self-similar, lattice- like or Bloch-like wave 
functions regarding the nature of the eigenstates. The associated spectra can be absolutely 
continuous, singular continuous or point singularities, and in many cases the spectra seem 
to be a mixture of these possibilities. All these information lead to a somewhat confusing 
situation, and these controversies can be resolved only if one has clear knowledge of every 
eigenvalue and eigenfunction of the system. The analytical answer to this problem is not 
very straightforward and in general, every quasiperiodic or aperiodic system needs separate 
attention. 

As a first step, in this paper we shall provide an instance of an aperiodic system where 
we have succeeded in overcoming all the ambiguities in its electronic properties. This is the 
case of the well-known Thue-Morse lattice (TM), a deterministic aperiodic system whose 
Schrodinger equation has been analytically solved to give a full account of all its eigenval- 



ues and eigenf unctions. Recent investigations on this system show that it exhibits many 
contrasting physical behaviors. The electronic properties of this lattice is more like those of 
periodic lattices than the quasiperiodic systems, while with regard to Fourier spectrum it 
has a degree of order intermediate between disorder and quasiperiodic lattices. The Fourier 
spectrum is singular continuous whereas the phonon spectrum is Cantor-set like W5 



Even if we focus our attention only on the electronic properties of the TM lattice, there are 
a number of papers in the literature showing many apparently contradicting features. Using 
the dynamical trace-map |jl6| approach Ryu et al. |jl3| have provided numerical evidence 



for the existence of critical and extended states in the TM lattice within the tight-binding 
model. From their results it also appears that lattice-like and chaotic wave functions are 
possible in the system. La Rocca [0 performed a multifractal analysis of the system and 



suggested that the spectrum contains absolutely continuous parts (extended states) as well 



as point singularities (localized states). Huang et al. [|T^ have calculated the mean resistance 
of the TM system and claimed that it is more localized than the Fibonacci chain. In an 
interesting analytical work, Chakrabrarti et al. []TU[ have provided the restrictive conditions 



on the energy eigenvalues for which the TM lattice can support extended electronic states. 
They have shown that a definite type of short range correlation among the atoms is respon- 
sible for the existence of extended states in the TM lattice, and it is quite different from the 
well-known dimer-type correlation But this work does not provide any insight about 
the entire spectrum giving full details of all the eigenstates of the system. Thus it seems 
that a clear understanding of the electronic properties of the TM lattice is still lacking and 
the present paper intends to resolve it on the basis of rigorous mathematical and physical 
considerations. 

The formalism we introduce is based on the dynamical trace-map [|T^ of the TM lat- 
tice and the commutation properties of the transfer matrices, and it provides a very simple 
analytical method for finding all the eigenvalues and eigenf unctions of the system. This is 
essentially a renormalization group approach and it gives a transparent picture of the nature 
of the states. The key point of this work is the association of the eigenvalues obtained by 



the trace-map technique with the commutation properties of the transfer matrices and the 
underlying physical mechanism is as follows. We can consider the infinite TM lattice, at any 
length scale, as built up of two basic clusters and the transfer matrices for these clusters be- 
come identical at certain eigenvalues of the system. So we can map the aperiodic TM lattice 
effectively into a periodic lattice, and it explains why this lattice can support extended elec- 
tronic eigenstates. Most surprisingly we observe that the TM lattice sustains only extended 
wavefunctions, a feature quite uncharacteristic of any quasiperiodic or aperiodic system. 
Also, we have shown exactly that the localization length ^ is always greater than the system 
size, which again signifies the delocalized character of all the states. Quite interestingly we 
have noticed that the symmetry of the TM lattice manifests itself in a beautiful way giving 
rise to the appearance of degenerate states in one-dimension within the tight-binding for- 
malism. This a very unusual behavior for any one- dimensional tight-binding system where 
degenerate eigenstates are almost forbidden, and the TM lattice can be considered as an 
exception. We have seen that apart from the four eigenvalues corresponding to the global 
band edges, all other eigenvalues are doubly degenerate, and thus our work invalidates the 



findings of Cheng et al. WM. To gain further insight about the nature of the states, we have 
also studied the transport properties of the electrons. We have derived exact expression 
for the Landauer resistivity of the lattice and shown that all the degenerate states exhibit 
Bloch-like behavior in the sense that the Landauer resistivity are zero for these states, while 
for the remaining four states the Landauer resistivity scales as ~ L^, L being the length of 
the chain. 

The TM sequence a"„(cr„) of order n may be generated by starting with the symbol A{B) 
and applying the inflation rule A AB and B BA n-times recursively. It is easy to 
see that o"„ is the complement of cr.„ obtained by interchanging A and B in cr„, and they 
satisfy the recursion relations (Xn+i = <^n(^n and o"„_|_i = o"„cr„. Now one can construct a 
TM chain by considering the symbols A and B representing two atoms A and B. The 
electronic properties of such a lattice are generally modelled by the on-site version of the 
nearest neighbor tight-binding Hamiltonian and the corresponding Schrodinger equation in 



Wannier basis is given by 

{E - ei)tjji = iJi+i + tjji_i , (1) 

where ipi is the amphtude of the wave function at the zth site, = or e^, and t is the 
hopping integral being set to unity. Now first we solve the Eqs. for the nth generation 
chain with the periodic boundary condition, and then consider the limit n — > oo for finding 
the electronic structure of the infinite lattice (see Ref. [|l9l). Using the transfer matrix 
method we now show that all the eigenvalues of the TM chain can be obtained from its 
trace- map which satisfies the relation 

= 4a^_i(a„ - 1) + 1 for n > 2, (2) 

where a„ = (TrM„)/2, M„ being the transfer matrix of the nth generation TM lattice. If 
we choose the origin of energy in such a way that = = V , then the initial conditions 
are ai = {E"^ — — 2)/2 and 0:2 = 2a1 — 2V'^ — 1. We see that a„ is a polynomial in E 
of degree 2", and for periodic boundary condition the relation an = ^ determines all the 2" 
eigenvalues of the nth generation chain consisting of = 2" atoms. Let us now recast the 
trace-map Eq. (Q) into the following form 

An-2 2 2 2 2/ i\ , i /'q\ 

This form of the trace-map is quite suggestive and gives a lot of information about the 
spectrum. The eigenvalues of the system are given by the roots of the equations 

al = al = ■ ■ ■ = a^_2 = and 0:2 — 1 = 0. (4) 



The equation 0^2 — 1 = actually determines the global band edges EO] , because the energies 



for which 02 > 1 are disallowed as it implies that a„ > 1. The global band edges are situated 



at E = ±(vKM-T ± 1), and Eqs. (^) show that these are the only four non-degenerate 
eigenvalues in the spectrum of the TM chain. It is clear from the rest of the equations in 
(I) that all the other eigenvalues are doubly degenerate. The form of Eq. (||) also suggests 



that these degenerate eigenvalues are all distinct. This is due to the fact that a„ cannot 
have any of the earlier a^'s (with i < n) as its factor. From the structure of Eq. it is 
obvious that the spectrum of the nth generation TM lattice contains all the eigenvalues of 
every preceding lower generation chains. In other words, the eigenvalues of any arbitrary 
generation TM lattice remain as eigenvalues in all the succeeding higher generation chains, 
and by the method of induction they also belong to the spectrum of the infinite TM lattice. 
In practice, the spectrum of the infinite lattice can be obtained by making the generation 
number n sufficiently large, and we get two global bands (symmetric around E = 0) where 
all the eigenvalues are doubly degenerate excepting those corresponding to the four global 
band edges. 

Knowing all the eigenvalues of the TM lattice, let us now look into the nature of the 
eigenstates of the system. From the symmetry of the TM sequence, we observe that the 
infinite lattice can be built up from every pair of unit cells of the form {an, ct„), where an{(Tn) 
is the nth generation TM sequence. If we define the corresponding pair of transfer matrices 
as (T„,T„), then in the basis of (2 x 2) identity matrix / and three Pauli matrices a^, Uy 
and az both T„ and T„ can be expressed into the following form: 
For odd n (n > 3), we have 

(5) 

where the coefficients satisfy the recursion relations 

«n = »n-2 + Pt-2 + ln-2 + + 2(7n-2 + ^1-2) 

X(3a^„2 + 1^1-2) + 2(7n-2'^n-2 " "n-2/^n-2) > 

I3n = 8a„_2/3„_2(7L2 + ^n-2) > (6) 

7n = 4(a^_2 - (3l_2 + 7n-2 + 6l-2)(^n-2ln-2 , 
= 4(a^_2 - + ln-2 + 5^_2)ttn-25n-2 , 

and the iteration starts from the set ai = 2ab — 1, /3i = 6 — a, 71 = —i{a + b) and 61 = 2ab, 
where a= {E-V)/2 and b= (E + V) /2. 
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For even n (n > 4), we get 

and the recursion relations of the coefficients are given by 

(«n-2 + li-2 + '^n-2)(«n-2 + 7l-2 + C-2) ) 
7„ = a„-2(6'n-2Atn-2 + 0n-2'^n-2) , 

(^n = «n-2(0n-2Ain-2 " ^n-2'^n-2) , (8) 
7^ = a„-2(6'n-2/Un-2 " 0n-2'^n-2) , 
5^ = a;„_2(0„-2/in-2 + Qn-2Vn-2) , 

where 6I„_2 = 2(7„_2 + 7^.2), K-2 = 2((5„_2 + S'^_2), /i„-2 = al-2 + ln-2in-2 + ^n-2'5L2 

and Vn-2 = ln-2^'n~2 ~ l'n~2^n-2- Here the initial set for iteration is 0^2 = — + 7^ + 5^, 
72 = 2(ai7i - 82 = 2{aiSi + i/9i7i), 7^ = 2(ai7i + i/^i^i) and 5'^ = 2{aiSi - ipi'ji). 

First we consider the case of odd n. If we set /5„ = in Eqs. (|^), then the matrices 
Tn and T„ become identical and this is exactly equivalent to the commutation relation 
[Tn-i,Tn-i] = 0. The condition /3„ = implies that Pn-2 = Pn-4 = ■ ■ ■ = = ^ {Pi 
being non-zero) and thus it also guarantees the equality between and for every 
m = 3, 5, 7, ■ ■ ■ n. It can be shown from Eqs. (|) and (j^) that the condition /?„ = 
actually satisfies the eigenvalue equation a„ = 1. In fact, we can reduce /?„ = as the set 
of equations ai = = = ■ ■ ■ = a„_2 = and 7i + = 0, and this is effectively a 
subset of the eigenvalue equation Eqs. obtained by the trace-map technique. It may be 
noted that the equation 7^ + 5^ = is identical to the relation 0:2 — 1 = 0, which gives the 
global band edges. Now we show that the rest of the equations in actually correspond 
to the equality of T„ and T„ for even values of n. From Eqs. (|^) and (|) we see that T„ 
and Tn become identical when a„_2 = 0, and this condition also satisfies the eigenvalue 
equation a„ = 1. Finally, combining the results for odd and even n, we obtain the equations 
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«i = «2 = ■ ■ ■ = «n-2 = and '-yf + Sf = as the conditions for the equahty between 
Tm and T^, where m = 3, 4, ■ ■ ■ n. It should be noticed that these equahty conditions 
precisely give all the distinct eigenvalues of the system (see Eqs. (^)), though they cannot 
detect the degeneracy of the levels. Nevertheless, it actually proves that at every eigenvalue 
of the system, we can always find a pair of unit cells for the TM lattice which become 
identical so far as the electronic properties are concerned. Thus the aperiodic TM lattice 
effectively maps into various periodic lattices for energy eigenvalues of the system, and this 
is the physical reason why all electronic states of the TM lattice are extended in nature. The 
wave functions do not have any Bloch-like periodicity, but they are extended in the sense 
that the amplitudes do not decay at infinity. 

For the sake of illustration, in Fig.l we display \ipn\'^ as a function of the site index n 
for a few selected energies. Here we have used the initial conditions ipi = ipQ = 1 which 



is consistent with the periodic boundary condition [|T^], and the figures clearly indicate the 
extended character of the wave functions. It is worth mentioning that we do not find any 
chaotic wave functions as reported in Ref. ||T3[. Normally the critical wave functions show 
power-law divergence [|T|, however, in the case of TM lattice all the states are well-behaved 
without any divergence at infinity. In Ref. |2T[] K. Iguchi has shown that the trace map Eq.(^ 
has no cycle, which implies that the TM lattice does not support power-law diverging states. 
It is also apparent from the above analysis that the transition from extended type states 
to critical or localized type states cannot be induced in the TM lattice by changing the 
Hamiltonian parameters, a behavior analogous to that of periodic systems. 

Now we calculate the Landauer resistivity p of the TM lattice. Since we can express 
the global transfer matrix M„ in the form (where M = unit cell transfer matrix and 
m = total number of unit cells) corresponding to every eigenvalue of the system, it can be 
written in the close form 
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Mn = I (for«„ = 0, n = 1,2,3, ■■■) , 

(for 72 + 51 = 0) 



r • \ 

1 + mos —tTwys 



y irwys 1 - m^s y 

By embedding a finite size TM chain in a periodic array of A type atoms connected by 
hopping integrals t = 1, we obtain the following expression for the Landauer resistivity [0, 

p = (when a„ = 0, n = 1, 2, 3, ■ ■ ■) 
= m^^l tan2(A;/2) (when 61 = i'ji) 
— m^S^cot'^{k/2) (when 61 = —i'ji) , 

where the wave vector k is given by the relation 2tcosk = E — e^. Thus we can further 
classify the states according to the transport properties of the electrons in the TM lattice. 
We see that all the degenerate states are Bloch-like in the sense that the Landauer resistivity 
is zero for these states, while the four states corresponding to the global band edges are not 
Bloch-like as in these cases p ~ L^, L being the length of the chain {L = N = ma, a = 
unit cell length). The non-zero value of p at the band edges physically indicates a transition 
between exponentially localized (gap) and extended (band) regions 0]. 

Let us now calculate the localization length (^) of the TM lattice from the expression 



for the global transfer matrix M„. The Lyapunov exponent can be evaluated exactly 
from the formula 7 = lim^y-^ooi^""'^ In ll^nll]? where || || is the modulus of the matrix. We 
observe that for all the states of system size A^, the localization length always has the value 



^ = 1/7 = N/ In 2, and it proves the absence of localization in the TM lattice [^. So from 



the point of view of localization of the electrons all states behave identically, whereas in 
terms of their transport properties the band edge states behave quite differently from the 
other states. 

To conclude, in this paper we have developed an analytical method which gives full 
information about the electronic structure of the TM lattice. We have shown that all the 
electronic states in the aperiodic TM chain are extended in nature, which is a very surprising 
and completely new result. Another interesting finding is that it is an instance of a tight- 
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binding one-dimensional system which supports degenerate eigenstates. In this formalism it 
is also possible to determine the Landauer resistivity and localization length of the system 
exactly. We hope that this method will provide a guideline for studying the electronic 
structure of systems of similar type. 
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FIGURES 

FIG. 1. Plot of IV'nP versus site number n. Here V = 0.5 and (a), (b) and (c) correspond, 
respectively, to energies 2.1180339, 1.5 and 1.9955076. All energies are measured in units of hopping 
integral t. 
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